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Abstract. Given a chain of viscoelastic spheres with fixed masses of the first and 
last particles. We raise the question: How to chose the masses of the other particles 
of the chain to assure maximal energy transfer? The results are compared with a 
chain of particles for which a constant coefficient of restitution is assumed. Our 
simple example shows that the assumption of viscoelastic particle properties has 
not only important consequences for very large systems (see |l]) but leads also to 
qualitative changes in small systems as compared with particles interacting via a 
constant restitution coefficient. 



1 Introduction 

We consider a linear chain of inelastically colliding particles of masses m^, 
radii Ri (i = . . .n), with initial velocities vq ~ v > and Vi = {i = 1 . . . n) 
at initial positions Xi > Xj for i > j with Xi^i — Xi > Ri+i+Ri. The masses of 
the first and last particles mo and m„ are given and we address the questions: 
How have the masses of the particles in between to be chosen to maximize the 
energy transfer, i.e., to maximize the after-coUisional velocity v'^^ of the last 
particle. If n is variable, how should n be chosen to maximize w^? Throughout 
this paper we assume that the initial distance of the particles is large enough 
to neglect "multiple collisions", i.e., only the first impact of each particle 
influences the final velocity of the iVth sphere of the chain. 

Recent investigations show that the properties of very large systems of 
viscoelastic particles differ significantly from those of particles interacting 
with constant coefficient of restitution . The system considered here may 
serve as an example of a small system which properties change qualitatively 
when the viscoelastic properties of the particles are considered. 

The coefficient of restitution is defined via e = | {v'^^i — v'^ / (f i+i — Ui)], 
which relates the relative velocity of the particles after the collision to the 
pre-coUisional quantity. The elastic collision corresponds to e = const. For 
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this case, basic mechanics yields: 

~^ mo , „ V 

The final velocity v',^ is maximized by 
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The optimal mass distribution for the case of a constant coefficient of 
restitution e does not depend on the value of e and is, therefore, the same as 
the optimal mass distribution in a chain of elastic particles. Figure ^ (left) 
shows the optimal mass distribution for different chain lengths n. The mass 
of the first particle is too = 1 and of the last particle to„ = 0.1. 




Fig. 1. Left: Optimal mass distribution rrii, i = 1 . . . n, for the case of constant e. 
Each of the lines shows the mass over the index i for a specified chain length 
n. The masses of the first and last particles are fixed at mo = 1 and m„ = 0.1. 
Right: Velocity distributions of the particles in chains of length n with the optimal 
mass distribution according to (|[) as a function of the chain length. The dissipative 
constant isti = l — e = 5- 10~*. The last particle reaches its maximal velocity for 
chain length n* = 44 (bold drawn). The velocity of the first particle of the chain is 
Vo = 1. 



In contrast to the mass distribution the corresponding velocity distribu- 
tions do depend on the value of the restitution coefficient e. Figure ^ (right) 
shows the velocity distribution for 6 = 1 — e = 5 • 10"'*. For the case of dissi- 
pative collisions the ratio — v'j-^/vq does not monotonously increase with 
n as for elastic particles (e = 1), but rather it has an extremum which shifts 
to smaller chain lengths with increasing dissipative parameter b. The optimal 
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value of n, which maximizes i?„ reads 

n* = log (mo/m„) / log (xq) (3) 
where Xo is the solution of the equation 

(l + xo) = (l + eK"/('+^°^ (4) 
Correspondingly, the extremal value of the Ry reads 

" 1 + e 



n„ — 



l+Xo 



(5) 



2 Chains of viscoelastic particles 



It has been shown that for colliding viscoelastic spheres the restitution coeffi- 
cient depends on the masses of the colliding particles and also on their relative 
velocity Vij n.An explicit expression for the coefficient of restitution is given 
by the seriesg,0 (see also 0) 
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with Y and v being the Young modulus and the Poisson ratio, respectively 
and = R,R,/{R, + R,), m^^ = m,m,/(m, + The material constant 
A describes the dissipative properties of the spheres (for details see ) . The 
constants Ci = 1.15344 and C2 = 0.79826 were obtained analytically in 
Ref. jj] and then confirmed by numerical simulations. 

In the following calculation we neglect terms O (w^/^) and higher and 
assume for simplicity that all particles are of the same radius i?, but have 
different masses. We abbreviate 



(8) 



Hence, for viscoelastic particles the velocities of the fc-l- 1-rst particle after 
colliding with the fc-th reads 
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Fig. 2. Optimal mass distribution of collision chains over the chain length n. The 
dissipative constant was 6 = 5- 10"^ (left) and b — 2 ■ 10""^ (right). 

The masses rrik , fc = 1 . . . ?i — 1 which maximize w'^ can be determined nu- 
merically and the results are shown in Fig. || for two different values of the 
dissipative constant b. 

For small chain length or small 6, respectively, the optimal mass distri- 
bution is very close to that for the elastic chain as shown in Fig. |l|. Again 
we find a monotonously decaying function for the masses. For larger chain 
length n or larger dissipation b, however, the mass distribution is a non- 
monotonous function. The according velocities of the particles in chains of 
spheres of optimal masses are drawn in Fig pi 




Fig. 3. The velocities of particles in optimal chains according to Fig. 
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2.1 Optimal mass-distribution 

The "loss" of energy, i.e., the amount of kinetic energy which is not trans- 
ferred from the first particle of the chain to the last one may be subdivided 
into "inertial" and "viscous" losses. Inertial losses occur due to mismatch of 
subsequent masses, which causes incomplete transfer of momentum even for 
elastic collisions if the masses differ. Viscous losses are caused by the dissipa- 
tive nature of collisions. The inertial loss is, thus, given by the energy which 
remains in the i — Irst particle after the collision with the ith: 
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We describe the chain in continuum approximation m(x) with rrii « Wi-i + 
'^"^^'^ ■ 1, where we assume that particles are separated on a line by unit 
distance. Discarding high-order mass gradients within the continuum picture 
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1 we write for the "line-density" of the inertial losses 
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Viscous losses may be quantified as the difference of the kinetic energy of 
a particle after an elastic collision and that of after a dissipative collision: 
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Expanding (|T^) up to linear order in the dissipative parameter b which is 
assumed to be small and neglecting products of b and mass gradients (which 
are supposed to be small too), the continuum transition of Eq. (12) yields 



dE„ 



,^3/5^11/5^ 



dx 23/5 

Thus, the total energy loss in the entire chain reads 
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with nix = dm/dx. For the second part of Eq. ( [T^ we assume in zero- 
order approximation the "ideal chain Ansatz" for the velocity distribution 
v{x), which refers to the velocity distribution v(x) in an idealized chain, 
where the kinetic energy completely transforms through the chain, i.e., where 
^m{x)v'^{x) = const — ^mow§. With rriQ = 1, vq = 1, i.e., v{x) = l/^/m(x), 
the right hand side of Eq. (jl^) follows. 

The mass distribution which minimizes Etot satisfies the Euler-equation 
applied to the integrand in (JT^): 



d 2mx d 
dx 8rn^ dm 
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(15) 



which leads to an equation for the mass distribution of the optimal chain, 
written for y{x) = l/m(x): 



dx^ 



(16) 



Figure 1^ (lines) shows the numerical solution of (p^, i.e., the optimal mass 
distribution, for a chain of length n = 40 for different damping parameters b. 
The points show the results of a discrete numerical optimization of the full 
chain problem (see Eq. (^)) applying a steepest descent method to optimize 
the masses rrifc of all particles. For small dissipation b both results agree. 




20 30 
particle number 



20 
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Fig. 4. Left: Mass distribution of the optimal chain of length n = 40 for different 
values of the dissipative parameter b. Lines: numerical solution of Eq. (|l^), Points: 
discrete numerical optimization (from top to bottom: • : b — 0.128, ■ : fe = 0.064, 
♦ : 6 = 0.032, A : 6 = 0.016, M: b ^ 0.008, T : 6 = 0.004, ►: 6 = 0.002, etc.). 
Right: Same data and symbols as left but plotted in larger scale. 



For larger values of b the solution of Eq. (|l^) deviates from the discrete 
optimization which is understandable since in our approximation we assumed 
the gradients of the mass distribution to be small which is violated for larger 
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b. While the absolute values of masses deviates from the discrete calculation, 
Eq. ( p^ still predicts well the position of the maximum of m^x). 

Figure |^ displays the corresponding distribution of velocities for the op- 
timal chains shown in Fig. ^. According to the maximum in the mass distri- 
bution, the velocity distribution reveals for larger b a pronounced minimum. 




10 20 30 40 1 2 30 4 5 60 

particle number nb^^" 

Fig. 5. Left: The velocity distribution along the optimal chain shown in Fig. ^ 
Lines from top to bottom: 6 = 2.5 ■ 10"", 5 • 10"", 0.001, 0.002, 0.004, 0.008, 0.016, 
0.032, 0.064, 0.128. 

Fig. 6. Right: The mass of the heaviest sphere m* in an optimal chain depends 
on the dissipative parameter b and on the chain length n. In the figure we plotted 
(m*)^^" over n\/6 for about 3000 different combinations of b and n (n — 2 . . . 300, 
b = 0.0001 . . . 0.256) including all data presented in Figs. ^ ^. Without any ad- 
justable parameters the data from the numerical optimization of chains agrees well 
with the analytical expression Eq. (p6|). 



One can give a simple physical explanation of the appearance of a maxi- 
mum in the mass distribution (and the corresponding minimum in the veloc- 
ity distribution) : As it is seen from Eq. ^ the restitution coefficient increases 
with decreasing impact velocity and increasing masses of colliding particles; 
this reduces the viscous losses. Thus, slowing down particles, by increasing 
their masses in the inner part of the chain, leads to decrease of the viscous 
losses of the energy transfer. The larger the masses in the middle and the 
smaller their velocities, the less energy is lost due to dissipation. On the 
other hand, since the masses ttiq and m„ are fixed, very large masses in the 
middle of the chain will cause large mass mismatch of the subsequent masses 
and, thus, large inertial losses [see Eq. (jl^)]. The optimal mass distribution, 
minimizing the total losses, compromises (dictated by b) between these two 
opposite tendencies. For the case of a constant coefficient of restitution the 
relative part of the kinetic energy, which is lost due to dissipation does not 
depend on the impact velocity. This means that only minimization of the 
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inertial losses, caused by mass gradient, may play a role in the optimization 
of the mass distribution. Thus, only a monotonous mass distribution with 
minimal mass gradients along the chain may be observed as an optimal one 
for the case of the constant restitution coefficient. 



2.2 The maximum of the optimal mass-distribution 

The mass m* of the heaviest sphere in the optimal mass distribution can be 
expressed as a function of the chain length n and the dissipative parameter h. 
With the term (dy/dx)'^ discarded, Eq. (^6|) describes formally scattering 
of a particle of unit mass by the potential 

U{y) = --dby^/^ with ^=-2^/^ (17) 

Formally changing notations x ^ t ("time") to emphasize the mechanical 
analogy, we write the equation of motion 

.. dU , , 1 / ^ 1 

y^--r^ 2/o = yi = = , y„^y{t^n)^ . (18 

dy TTiQ rUn 

Here we consider the case of mass distributions having a maximum; the gen- 
eralization, however, is straightforward. Hence, 

+ U{y) = const - U[y*) (19) 

where y* is the turning point in the scattering problem, i.e., the point where 
the particle's "velocity" y is zero (this corresponds to the point m* of the 
mass distribution in the initial problem). The "particle" reaches this point 
at "time" t* , i.e., 

f^dh(y'l^-y*^'^) . (20) 
Solving Eq. ( po|) with respect to y yields 

|=±Vd6y*^/V(y/r)'''-l- (21) 

Integration over "time" from t = to i = n in Eq. (^l|), therefore, leads to 
(with correct choice of signs) 
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Using the substitute z = (y* lyf^'^^ the integrals in Eq. ( p2| ) may be recasted 
into the form 



z 10 (1 — z) dz = B 
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(23) 



with fc = for the first integral in the the left-hand side of Eq. (g^) and with 
k = n for the second integral. B{x,y) is the Beta-function and B{x,y, a) is 
the incomplete Beta-function (which has an upper limit o instead of 1 in 
its integral representation). If we assume the pronounced maximum in the 
optimal mass-distribution, so that a = {y* /ykf'^'^ = (mfe/m*)^^^ is small, 
one can approximate the incomplete Beta-function as 



B I — , -, a 

10' 2' 



z 1" {l — z) ^ dz 



lo Jo 
With the use of Eqs. (||) and (p|), Eq. (|2|) reads 
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In the original variables, which refer to the mass distribution, we obtain a 
scaling relation, connecting the heaviest mass m* , the chain length n and the 
dissipative parameter b: 
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with the constants 
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So far we considered the solution of the variational Eq. (|l^) with the term 
y^^ (dy/dx) discarded. For this case the constant d, which has been given 
above reads: = | 2^/^. It may be shown, however, that perturbative (thus 
approximate) account of this omitted term leads to an equation of the same 
form as Eq. (pT]), but with the renormalized coefficient d | d = J|2^/^; 
the details are given in Using numerical values for B (^j^,^) (see the 
renormalized coefficient d and mg — 1, m„ = 0.1 for the first and the last 
masses of the chain, yields for p and q: 



p = 0.15217 q = 0.68989 (28) 

In Fig. ^ we compare the analytical relation Eq. ( |2^ ) with the constants 
given in Eq. (E^) with numerical results for to*. The numerical data follow 
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from the numerical optimization of the mass distribution for different chain 
length and different dissipative constants, including all data given in Figs. |^ 
||. As one can see from Fig.^, the results of the analytical theory and of the 
numerical optimization agree well, except for large dissipation values. We 
would like to stress that no fitting parameters have been used. 

3 Conclusion 

We investigated analytically and numerically the transmission of kinetic en- 
ergy through one-dimensional chains of inelastically colliding spheres. For 
constant restitution coefficient, e = const., the distribution of the masses 
which leads to optimal energy transfer, is an exponentially decreasing func- 
tion which is independent on e, i.e., it is the same as for elastic particles with 
e = 1. 

For viscoelastic particles where e depends on the impact velocity, the 
optimal mass distribution is not necessarily a monotonous function, but de- 
pending on the chain length n and on the material parameters of the spheres 
it may reveal a pronounced maximum. 

We develop a theory which describes the total energy losses along the 
chain, so that the optimal mass distribution, minimizing the losses, may be 
obtained as a solution of a variational equation. We derived an expression 
relating the heaviest mass in the chain to the chain length and the dissipation 
constant. Having no fitting parameters, it is in good agreement with the 
numerical data. 

It has been demonstrated before that for the case of "thermodynamically- 
large" granular systems the impact-velocity dependence of the restitution 
coefficient leads to qualitatively different behavior as compared to systems 
with e = const, (e.g. |p]-p|,|lO|-[l^). Our system demonstrates that the veloc- 
ity dependence of the restitution coefficient leads to qualitative modifications 
in small and simple systems too. Therefore, in general, we believe that the 
assumption of a constant coefficient of restitution is an approximation which 
justification cannot be assumed a priori but has to be checked for each par- 
ticular application. 
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